The Fastest Rotating Pulsar: A Strange Star?* 
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According to the observational limits on the radius and mass, the fastest 
rotating pulsar (PSR 1937+21) is probably a strange star, or at least some 
neutron star equations of state should be ruled out, if we suggest that a dipole 
magnetic field is relevant to its radio emission. We presume that the millisecond 
(*C) ! pulsar is a strange star with much low mass, small radius, and weak magnetic 

moment. 



PACS: 97.60. G, 97.60.J, 95.30.C 



Pulsars are conventionally modeled as neutron stars, but a strange star0 i model 
for pulsars is also proposed (see also a review, e.g., by Xu et al.il). Therefore it is of 
great importance to find observationally unambiguous features to distinguish strange stars 
from neutron stars. There are three hopeful ways known hitherto to identify a strange 
star, which are based upon the differences of viscosity, mass-radius relation, and surface 
condition between neutron and strange stars, cl Recently, Kapoor and Shukrei suggested 
to constrain the equations of state of neutron stars by Rankin's experiential linel of the 
core emission of radio pulsars, and find the restriction that the pulsar masses M < 2.5M Q 
and radii R < 10.5 km, indicating that puslars are strange stars. However, Kapoor and 
Shukre's wor kl has at least two unseemly points: (1) the polar cap is defined for aligned 
rotators (rather than for orthogonal rotators) in their calculation; (2) the Rankin line is 
not doubtless true due to many observational and statistical uncertainties. By redefining 
the polar cap for orthogonal rotators, still we can not obtain a conclusion that puslars are 
strange stars, or that some neutron stars equations of state are ruled out, if the Rankin line 
is used to constrain pulsar masses and radii. Nevertheless, when applying this method to 
the fastest rotating pulsar, PSR 1937+21, we conclude that the pulsar may be a strange 
star, or at least some equations of state should be ruled out. In order to find solid evidence 
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to shown that PSR 1937+21 is a strange star, it is strongly suggested to measure precisely 
the maximum rate of position-angle swing by further polarization observations. 

Now, we improve Kapoor and Shukre's computation by removing the first unseemliness. 
In a spherical coordinate system with magnetic axis /i being chosen as the z-axis, one has 
the following form for a position vector r in a dipole field line denoted by parameter A,@ 

r = A— sin 2 #, (1) 
2n 

where 9 is the polar angle, P the rotation period, c the speed of light, and A characterizes 
the sorts of field lines. Based on Eq. ([]]), the angle 6* M between \i and the direction of 
magnetic field at r reads, 

2A-3^ 

cos 6> M = —= (2) 
'A(4A - m) 



where 1Z = ^ <C 1 for points near pulsar surface. Expand function 6^ in the vicinity of 
1Z = 0, one comes to 

9 ^^ I/2 + ^ 3/2 + 0(K6,2) - (3) 

Assuming A = Ao for the last-open-field lines (note: Ao is a fuction of inclination angle a), 
one can obtain the beam radius p for the radiation at an emission height h = r — R (R is 
the radius of pulsar), neglecting the terms being equal or higher than 1Z 3 ^ 2 in Eq.(|3|), 

1 24° 

" ' (4) 
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For core emissions (r ^ R), which are supposed to be originated from regions 
near pulsar surfaces, general relativistic effects are not negligible due to the spacetime 
curvature.'! Two such effects, squeezing of the dipole magnetic field and bending of the 
radio wave, can be represented approximately by two factors, S i / sqz and /bnd, respectively, 
so that the beam radius p can be re-written as 



_ 3 /2tT r 1 ' 2 L24° . , 

P — TPTF? ~ A— r>\/ TTvi isqz/bnd, \?) 



2 V XocP 1 ' 2 10km JsqzJbnd ' 

where / sqz = (1 + f^)~ 1/2 , And = |(2 + (1 - ^f)^ 1/2 ), G is the gravitation constant, M 
the pulsar mass. 

We are to calculate the parameter Ao of last-open-field lines below, assuming that only 
the plasma within the light cylinder with radius r\ c = cPj (2ir) can corotates with the star, 
i.e., field lines that penetrate beyond the light cylinder should be open. Generally, Ao is a 
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function of inclination angle a. Nevertheless, one can obtain simple formulae of Ao in two 
special cases. Obviously, A (a = 0) = 1 for an aligned rotator according to Eq. ([[]). Another 
formula of Ao can be easily reached for an orthogonal rotator (a = 90°). The conal angle 8q 
of the null surface, where the magnetic fields are perpendicular to magnetic axis p, can be 
obtained from Eq.(§) by setting 6> M = 90°, 




Let rcos^o = r lc , one finds 

Ao(« = 90°) = ^5, (7) 

which is 2.6 times of Ao(ct = 0°). 
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Fig. 1. — Calculated results based on Eq . (|TTD show the variations of p x \fP (in deg-s 1 / 2 ) as 
a function of emission altitude r for different pulsar masses ranging from O.1M to ^ 2M . 
We present results of Ao = 1 (aligned rotators) and A = 3\/3/2 (orthogonal rotators) for 
two limits; an actual value for a pulsar with artificial inclination angle a should be between 
those two corresponding values. The dashed lines are calculated by Eq.([|), which are not 
relevant to pulsar mass. 

Following the representation by Kapoor and Shukre,!! we calculate the value, p x \f~P ', 
as a function of emission altitude r for pulsars with core emission according to Eq. ([□]). 
The calculated results are shown in Fig .[I]. Contrary to the result of Kapoor and Shukre,@ 
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it is clear from Fig.[l] that the Rankin line (i.e., p x \J~P = 1.225° s 1 / 2 ) can not effectively 
put any limit on the equations of state of neutron stars or conduce toward any indication 
that pulsars are strange stars if we assume that Rankin's samples are of nearly orthogonal 
rotators. Nevertheless, we may still obtain some valuable conclusions from above calculation 
for some special pulsars, for instance, the fastest rotating millisecond pulsar, PSR 1937+21. 

The rotation period of PSR 1937+21 is P = 1.558 ms, the smallest one of pulsars 
observed, and the rate of period change is P = 10 -19 s • s _1 . The radius of light cylinder 
of the star is r\ c = cPj (2tt) = 74km . The corotation radius, defined by the balance of 
gravitational force and the centrifugal force, is r c = [GM/(47r 2 )] 1 / 3 P 2 / 3 = 2O(M/M ) 1 / 3 km, 
where M is the mass of PSR 1937+21. Therefore the pulsar radius R satisfies 

R < min(r lc , r c ) = 2O(M/M ) 1/3 km. (8) 

Another circumscription of the stellar radius R and the mass M may arise from the 
inclusion of general relativistic effect. @ The pulse width A<p of PSR 1937+21 is about 10° 
at 600 MHz and 1.4 GHz (EPN database), which is difficult to be explained geometrically 
if we suggest PSR 1937+21 has a radius of 10 km for canonical pulsars. @ Gil propose 
thus a mo in which a quadrupole magnetic field geometry is assumed. However an 

alternative conjecture to overcome the difficulty is that PSR 1937+21 may have unusual 
small radius. Taking the simplest proposing that a = 90° and impact angle (3 = 0, we know 
Acj) = 2p, and therefore pVP = 0.2° s 1//2 . According to the results in Fig.|l|, it is found that 

M < 0.2 M© and R < 1 km. (9) 

These stringent limits have to result in the conclusion that PSR 1937+21 is a strange star 
rather than a neutron star because of the strikingly different mass-radius relations. 

However, the conclusion in Eq.([5J) is not so solid. In fact there are two models proposed 
to account for the inter-pulse emission: @ single-pole model and double-pole model. Three 
observational fact s§ favour a double-pole model of inter-pulse origin of PSR 1937+21. 
(1) The gradient of the position angle at inter-pulse center has the same value as at the 
main-pulse center. (2) The separation of the main-pulse and the inter-pulse is nearly 180°. 
(3) The intensities of the main-pulse and the inter-pulse are roughly equal. The observed 
maximum rate of position angle swing at the main-pulse center is {dip/d(p) maiX ~ 3. In the 
standard rotating vector model, the impact angle (3 should be 

/3 = sin- 1 [ S ""* ]~19.5°, (10) 
(d^/d0) max 

and the beam radius p thus can be derived as 

p = cos" 1 [cos/5-2sinasin(a + /3)sin 2 ^] ~ 20.1°, (11) 
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if we take a = 90°, (d^/d0) max = 3 and A0 = 10°. In this case, we find p\fP = 0.79° s 1 / 2 , 
and the limits thus are 

M<2.4M and R < 11.5 km (12) 

based on Fig.[l| 

It is well known that the masses and radii of rotating magnetized neutron stars are 
of equation-of-state dependent. The radius of a neutron star with a maximum mass is the 
smallest; a smaller mass, a slower rotation and/or a lower inner magnetic field would cause 
its radius to be larger. H3~ilJ Therefore the observational deduction of Eq.(|l2|) should have 



some implications for proof-testing the equations of state available. Five euaqtions of state 
have been focuse d@ in the calculation for rapidly rotating neutron stars, the results of 
which are listed in Table 1 by interpolating between the tabulated points in Table 9-23 of 
Cook et al.@ for angular frequency Q = 4 x 10 3 s -1 (i.e., P = 1.558 ms). We see that the 
equations labeled as "L" and "M" should be ruled out by the limits of Eq. (p~2"D . Certainly, 
we should also keep in mind that more equations of state may be killed by inclusion of the 
inner magnetic field effect, which is an interesting topic in the future study. 

Furthermore, two observational uncertainties may result in a derivation of stronger 
limits on the mass and radius of PSR 1937+21. The first is that observation would give a 
less steep position angle gradient (i.e., a smaller value of (d , 0/d0) max than reality) due to 
smearing of finite sampling time, to the frequency dispersion in pulse arrive time,U and to 
the interstellar scattering^ which may bring on a stretching out of the longitude scale. In 
fact, (d , 0/d0) max > 3. A larger (d-?/>/d(/>) max (thus a smaller (3 according to Eq.([LC])) favours 
a rough equality of the intensities of the main-pulse and inter-pulse. The second is that 
PSR 1937+21 might not be exactly an orthogonal rotator (i.e., a is not precisely 90°). It 
is possible that a < 90°. Both these two uncertainties lead to a smaller value of py/P. 
We thus take (d^/d0) max and a as two free parameters to calculate p\fP. The calculated 
results are shown in Fig|| It is clear that stronger limits than those of Eq . (|T2|) may rule out 
more neutron star equations of state, and probably conduce toward a strange star model 
for PSR 1937+21, based on Figs.|I] and @. 
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Fig. 2. — Beam angular radius times square root of period {p\fP) of PSR 1937+21 versus 
possible different inclination angles a for the case of (d-?/>/d0) max = 3, or maximum rates of 
position-angle swing (d?/>/d0) max for a = 90°. 

In Fig.|| we find p\f~P is sensitively dependent on (dvp/d(f)) maiX if (d^/d0) max ^ 10. 
p\[P = 0.6 if (d?/>/d0) max = 4 and a = 90°, in this case the limits of Eq.(|12|) should 
be modified as M < 1.4 M and R < 6.6 km. These limits inevitably lead to the 
conclusion that PSR 1937+21 is a strange star. Because of the observational uncertainties 
of (d^/d0) max and a and of the importance of distinguishing a strange star in nature, we 
seriously table a proposal of detecting more accurate polarization signal from PSR 1937+21 
to obtain an actual value, especially of (d^/d0) max . This could be possibly by means of 
reducing the sampling time and increasing the total observation time. 

If PSR 1937+21 is a strange star with small radiu (R < 10 km), the old formulae0 to 
calculate the surface magnetic field B and dipole magnetic moment p m should be modified 

as 

B = (^) 1/2 ~ 9.3 x W*R- ll \PP)^ G, 

^ m = {±£^t) V2 „ 4.6 x 10™ Rl /2 (Ppy/2 G-cm 3 , 1 ' 

since the density p ~ 5 x 10 15 g cm -3 has a very modest variation with radial distance of 
strange star, M ~ (4/3)7ri? 3 p, and the moment of inertia I ~ (8/15)7rpi? 5 , where R$ 
is the stellar radius in 10 5 cm (1 km). We see from Eq.(|IIJ) that the calculated magnetic 
moment is strongly related to the radius although the observation-determined surface 
magnetic field is weakly related to. For instance, if the radius of PSR 1937+21 is only 1 
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km, B = 1.2 x 10 9 G, whilst /i m = 5.8 x 10 18 G-cm 3 (typically, /i m for millisecond pulsars 
is assumed to be 10 26 G-cm 3 ). This property favors the assumption that dipole magnetic 
structure dominates in the radio emission region. 

Millisecond pulsars are currently believed to be of recycled-origin of normal pulsars 
which are spin down enough. However, motivated by the study of planet formation 
around PSR 1257+12, Miller and Hamiltonlll suggested that some and perhaps all isolated 
millisecond pulsars may have been born with high spin rates and low magnetic fields instead 
of having been recycled by accretion. This is understandable if we assume part or all of 
the isolated millisecond pulsars are strange stars with smaller radii (thus smaller masses): 
(1) Turbulent convection in nuclear matter!! or in strange quark matter© should be less 
prosperous in a proto-pulsar with lower mass than with higher mass in the Kelvin-Hermholtz 
cooling phase, the dynamo-created magnetic field is thus weaker. (2) Small radius favors 
faster rotation since the centrifugal force become smaller while the gravitational force to be 
larger (see Eq.flD). A bimodal distribution of pulsar periods and magnetic fields may arises 
from various kind of progenitors as well as complex mechanisms of supernova exploration. 
Therefore, it is conjectured that PSR 1937+21 and PSR 1257+12 (and possibly some or all 
of the isolated millisecond pulsars) may have weaker dynamo-originated magnetic field due 
to a less effective field magnification process. 
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Table 1. Total gravitational mass-energy M and the circumferential 
Radii R at equator of rotating (P = 1.558 ms) neutron stars with 
maximum mass for five kinds of Equations of state 



* See Table 2 of Cook et al.il for details. 



EOS* 


M/Mq 


R(km) 


A (Reid soft core) 


1.6604 


8.82 


L (Mean field) 


2.7263 


14.98 


M (Tensor interaction) 


1.8298 


19.08 


AU (AV14 + UVII) 


2.1433 


9.78 


FPS (UV14 + TNI) 


1.8069 


9.88 



